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FOREWORD 


Single  particle  motion  in  a  planar  undulator  wiggler  free 
electron  laser  magnetic  field  configuration  is  set  up  and  examined 
for  the  construction  of  invariants,  preparatory  to  specifying  beam 
equilibria  for  input  to  a  Maxwell-Vlasov  analysis.  Two  Invariants 
are  immediately  available,  the  energy  y  and  one  component  of  trans¬ 
verse  canonical  momentum  Py .  In  an  effort  to  find  a  third  invariant, 
or  a  basis  of  determination  for  an  advantageous  approximate  third 
invariant,  which  can  be  associated  with  the  axial  periodicity  of  the 
field,  a  variable  change  having  some  interesting  features  is  reported 
here.  The  axial  motion  between  z  =  +  ®  is  converted  to  the 
oscillations  of  a  particle  in  an  infinitely  long,  axially  reflecting 
box.  The  transit  of  a  half-cycle  of  wiggler  oscillation, 

|  Az  |  =  Tr/ko,  corresponds  to  a  single  bounce  between  Z  =  +  ®,  where  Z 
is  the  new  coordinate.  Transverse  coordinates  and  velocities  are 
continuous  at  each  axial  reflection  from  infinity. 
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CHAPTER  I 
INTRODUCTION 


Determination  of  beam  equilibrium  configurations  forms  the 
starting  point  in  analysis  of  the  free  electron  laser  instability 
for  an  intense  beam.  For  equilibrium  distributions  of  beam 
electrons  to  be  constructed  which  will  take  into  account  non- 
linearities  associated  with  beam  electron  motions  in  the  prescribed 
magnetic  field  through  which  the  beam  passes,  as  well  as  effects  of 
finite  transverse  geometry,  the  invariants  of  single  particle  motions 
in  the  given  field  are  needed.  When  beam  self-fields  can  be 
neglected,  the  determination  of  such  invariants  will  constitute  a 
solution  for  the  single-particle  motion;  the  equilibrium  distribution 
will  be  a  function  of  these  invariants.  Corresponding  modulational 
influences  on  beam  envelope  can  then  be  included  automatically  in  the 
stability  analysis. 

In  two  different  model  field  examples  this  has  been  done  by 
Bernstein  and  Hirshf ield, 1  Davidson  and  Uhm2>3  and  Uhm  and 
Davidson.^  In  Refs.  1-3  the  helical  wlggler  field  was  allowed  to 
have  only  axial  dependence. 


B 


wigg 


B  coskQz  x  -  B  sinkQz  y 


(1.1) 


where  kg  =  2ti/\q,  with  Aq  the  wiggler  wavelength,  and  where  B  is  the 

constant  amplitude  of  the  wiggler  field.  In  Ref.  3  a  uniform  axial 
guide  field  also  was  present.  In  Ref.  4  an  axial  undulator  wiggler 
with  a  guide  field  was  studied  for  the  case  of  a  hollow  beam,  which 
allowed  radial  variations  of  the  field  components  again  to  be 
neglected, 

B  =  Br  r  +  Bz  z,  (1.2) 

■4-  A 

where  the  axisymmetric  components  of  B  are  got  from  the  6  -  component 
of  a  vector  potential  given  by  A  =  A.  0, 
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Br  =  -  3A0/3z,  Bz  =  r-13(rA0)3r,  (1.3) 

with  A0  =  A0  (r,z),  and  with  also  v  x  B  *  0.  Setting  r  *  Rq,  the 

equilibrium  beam  mean  radius,  allowed  B  and  B  to  be  approximated  as 

r  z 

independent  of  r,  and  simplified  the  analysis  of  the  single  particle 
equations  of  motion.  The  field  then  has  a  form  similar  to  eq .  (1.1), 

Br  ~  "  ®r  sinkoz>  Bz  ~  B0  ”  Bz  C0skoz  (l-1*) 

where  Bq  is  the  uniform  guide  field,  and  Br,  B  are  constant  wiggler 
amplitudes  determined  by  Bq,  kgR  ,  and  the  mirror  ratio 
R  =  Bz  (r  =  0,  kQz  =  +  7r)/Bz(r  =  0,  kQz  =  0). 

In  Refs.  1  and  2,  three  exact  invariants  of  electron  motion,  the 
energy  and  two  components  of  transverse  canonical  momentum,  were 
available;  and  the  same  is  true  of  Ref.  3,  but  where  the  invariants 
(CA,  Ch,  C  ),  now  also  allowed  incorporation  of  helical  and  axial 

field  effects  together.  In  the  undulator  configuration  of  Ref.  4 , 
only  two  exact  Invariants  were  available,  the  energy  and  canonical 
axial  angular  momentum,  and  a  third  approximate  invariant,  the  axial 
component  of  particle  momentum,  pz,  was  employed.  So  the  undulator 
wiggler  field  single  particle  motions  have  also  proved  more  difficult. 

In  the  present  report  I  describe  a  brief  investigation  of  single 
particle  motion  in  a  planar  undulator  field,  specified  by  a  vector 

potential  A  =  Ay  (x,z)y.  Only  two  invariants  are  immediately 
available,  the  energy  and  y-component  of  canonical  momentum.  As  a 
starting  point  of  the  search  for  a  third  invariant  with  origin  in  the 
axial  periodicity  of  the  wiggler,  or  to  produce  a  convenient  basis 
of  determination  of  an  approximate  third  invariant  better  than  pz,  I 
have  transformed  the  axial  coordinate  z  to  a  new  coordinate  Z,  which 
ranges  between  +  °°  as  z  covers  a  half-cycle  of  wiggler  oscillation. 

The  transformation  is  fixed  by  specifying  that  the  wiggler 
amplitude  be  absent  from  the  axial  equation  (for  zero  guide  field). 
After  elimination  of  the  y  coordinate,  the  xZ-equations  describe  a 
particle  oscillating  in  a  box  of  infinite  "length"  and  having 
axially  reflecting  walls,  while  also  executing  its  transverse  motions. 
In  this  way,  the  full  motion  along  the  beam  axis  is  transformed  into 
motion  over  only  a  single  cycle  of  wiggler  oscillation. 

It  is  not  clear  whether  the  altered  problem  presents  special 
calculational  advantage,  and  time  has  not  permitted  more  than  a 
cursory  examination  of  this.  It  has  seemed  appropriate  nevertheless 
to  record  the  present  formulation  as  a  slightly  unusual  and  possibly 
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useful  starting  point  of  analysis  In  which  features  associated  with 
wiggler  periodicity  receive  explicit  attention  at  the  outset.  The 
method  is  more  generally  applicable,  e.g.  to  problems  of  cylindrical 
symmetry  as  well  as  to  the  present  case  of  planar  symmetry. 
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CHAPTER  II 

PLANAR  UNDULATOR  WIGGLER  EQUATIONS  OF  MOTION 


The  vector  potential  is  A  =  A  y,  and  the  magnetic  field's 
nonvanishing  components  are  y 


B 


x 


9A  3A 

_ L  B  =  —X- 

3z  *  Z  3x  * 


(2.1) 


The  choice 


Ay(x,z)  =  BQ  x  -  k"1  B1  sinh  kQx  sinkQz  (2.2) 

models  a  static  undulator  configuration  with  mean  guide  field  Bn 
and  undulator  amplitude  B, ,  where  Bn  and  B,  are  constants.  DueU 
to  the  symmetry  1  u  l 


Ay(-x,z)  =  -  Ay ( x , z )  (2.3) 

Bz  has  the  sign  of  x,  corresponding^to  like-pole  magnet  pairings 
successively  alternating  along  the  z  -  direction.  The  field 
satisfies  the  condition  V  x  §  =  0.  The  Hamiltonian, 

H  =  G  [Px  +  Pz  +  (Py  -  qc_1Ay(x,z))2  +  m2c2]  (2.4) 

where  c  is  the  speed  of  light  in  vacuo ,  q  the  charge  on  the  electron 
and  m  its  mass,  is  invariant  to  translations  along  the  y-axis  so  P 
is  a  constant  of  the  motion  along  with  H;  and  the  problem  of  y 

configuration  motion  is  therewith  reduced  to  a  two-dimensional  one 
in  x  and  z . 

It  will  be  enough  to  study  the  configuration  equations  of 
motion.  These  are 
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‘x-x>  i 


y '-^Bx  <x>*>  i 


(Vx-z)  x  -  Bz  (x-x)  x) 

■  -  3t  (  ymc  VX,Z>)  •  <2 

2 

the  last  by  the  constancy  of  y,  whose  value  is  H/mc  ,  and  by 
eqs.  (2.1).  Prom  eq .  (2.7) 

.  Py  qAyu,z)  (2 

^  ym  ymc  1 

where  Py  is  the  constant  of  integration;  whence,  upon  substitution 
into  eqs.  (2.5)  and  (2.6)  one  has 


d2x  _  _d£  JTz  _  d± 
dt2  '  '  sx  ’  at2  '  '  92 


where 


qPv  1  (q\)2 

*  -  *  (x,z)  =  -  -jV  A  + 

ymc  ' 


Eqs.  (2.9)  have  the  first  integral 


jr  (x2  +  z2)  +  i<i  (x,z)  *  C  =  const 


P  ,  y  and  C  satisfy  a  relationship 
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-yy 


(2.12) 


so  only  two  of  the  three  constants  are  Independent,  and  y  can  be 
exchanged  for  C  if  we  wish. 

Note  from  eq.  (.2.8),  P  =  0  corresponds  to<y>=  o  for  electron 

«7 

distributions  even  under  x  -*■  -  x  owing  to  eq .  (2.3)-  In  this  case 

p 

also,  2C  has  the  value  vn,  where  vn  is  the  constant  speed  of  the 
charge .  u  u 

I  assume  Bq  =  0  for  simplicity  now  and  expand  eqs .  (2.9)  using 
(2.2)  and  (2 


e  bo 
.10), 


d2x 


dt 


2  ~  ~  cosh  ICqX  •  sink^z  •  (a  +  b  sinh  kQx*sinkgZ) 


(2.13) 


d2Z 

— ~  *  -  sinh  knx *cosknz • (a  +  b  sinh  knx-sinknz), 
dt  u  u  v  u 


(2.1*0 


where 


a  =  qB1/(ym)2c  •  Py  =  ^  Fy 


(2.15) 


b  =  (qB1/ymc)5K0'1  =  k"1  p2. 


(2.16) 


with  =  qB-^/ymc;  from  eq.  (2.11)  also 


.2  .2 


-1 


x  +  z  +  2  a  kg  sinh  kgX^sink^z 


+  bk"1  sinh2  kgX*sin2kgZ  =  2C 


(2.17) 
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CHAPTER  III 
CHANGE  OP  VARIABLES 


If  a  new  axial  variable  Z  is  introduced. 


Z  =  k"1  g(kQz) , 


(3-D 


independent  of  t  so  the  Hamiltonian  will  be  unchanged  in  value,  and 
with  g  a  function  at  our  disposal,  we  have  the  Identities, 


2  g'(kQz) 


-1 


and 


(3-2) 


Z' (k0z)_1 


kQg"(knz) 


(3-3) 


where  primes  denote  differentiation  with  respect  to  the  argument. 

Eq .  (2.14)  becomes 

r  *2  -i 

Z  -  kgg  2C  -  x  -  2akQ  sinh  kQx.sinkgZ 
-12  2  1 

-bkg  sinh  kgX*sin  k^z  =  -  g'»a  sinh  kgX*coskQz 
2 

-  g'*b  sinh  kgX-slnkgZ  coskQz.  (3.4) 

Note  the  choice  g(fc)  =  6,  with  3  =  k_z  gives  Z  ■  z,  and  (3*4) 
reduces  to  (2.14)  as  it  should.  u 

The  x-dependence  of  the  coefficients  of  a,  to  be  called  a-terms, 
is  the  same  on  both  sides  of  eq.  (3-4);  the  same  is  true  of  the 
coefficients  of  b,  viz.  the  b-terms.  So  the  a-terms  can  be  trans¬ 
formed  away  by  a  suitable  choice  of  g;  or  the  b-terms  can  be 

13 
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removed  Instead.  (They  cannot  both  be  removed  by  a  transformation 
of  z.)  The  former  is  accomplished  if  g  satisfies 


2g"(e)  Sine  =  -  g'(  s)  cose  (3-5) 


i.e. 


g'  ( S) 


|  sine | , 


A  >  0. 


(3-6) 


The  transformation  cannot  be  determined  as  an  elementary  function 
owing  to  the  square  root  in  eq.  (3*6),  whose  cause  was  the  factor  2 
in  eq.  (3-5). 

For  the  b-terms  matters  are  different;  the  analogue  of  eq.  (3-5) 


is 


g"  ( B)  sin^e  -  -  g'(  b)  sine  cose, 

whence 

|g'(e)  sine |  =  | A | , 


(3.7) 

(3.8) 


or,  more  serendipitously , 

g' ( 6)  =  A  cose,  (3-9) 

where  A  is  an  arbitrary  nonzero  constant.  Integrating, 

kQZ  =  g(e)  =  A  log  |b  tan  ^  6],  B  =  const.,  (3-10) 


where,  it  will  be  recalled,  e  =  kgz.  B  may  be  assumed  positive 
without  loss  of  generality;  I  assume  further  that  B  =  1  as  there 
seems  to  be  no  gain  resulting  from  its  retention  as  a  free  parameter. 
(Figure  1  shows  a  sketch  of  g(fO  for  A>0,  and  0  <  f$  <  n . )  Solving 
eq.  (3.10)  gives  also, 

sine  =  +  sech  A'^qZ,  .  tane  *  +  csch  A'^k^Z,  (.3-11) 

the  sign  associations  for  fixed  A  depending  upon  the  angular  range  in 
which  B  falls.  From  (3. 11),  also, 
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cos0  *  -  tanh  A_1kgZ.  (3.12) 

Hence,  from  eqs .  (3.7),  (.3*9),  (3.11)  and  (3.12), 

g"(fj)  =  A  sinh  A”^kgZ*cosh  A  ^kgZ,  (3-13) 

which,  from  Inspection  of  eq.  (3.*0,  suggests  either  A  =  2  or 
A  ■  -  2.  Each  choice  gives 

g"(  $)=  +  sinh  kQZ,  |A|  =  2,  (3-12*) 

and  eq.  (3*^)  becomes 

Z  -  ( 2C  -  x2)  kg  sinh  kgZ 

*  +  2a  sinh  kQx*sinh  kgZ;  (3.15) 

from  eq.  (2.13) 

"1  2  1 
x  +  j  b  sech  -x;  kgZ*sinh  2kgX 

=  +  a  sech  h  kgZ*cosh  kgX.  (3.16) 

The  "energy  equation",  eq.  (2.17),  reads 

x2  +  z2  +  2akg!  sinh  kgX*sech  ^  kgZ 

+bk~^  sinh2kgX  •  sech2  ^  kgZ  =  2C,  (3.17) 

where 

z2  =  ^  sech2  ^  kgZ-Z2.  (3-18) 

The  upper  and  lower  signs  in  eqs.  (3*15)  -  (3-17)  correspond  to 
those  in  eqs.  (3*11).  Notice,  finally,  the  periodic  functions  of  z 
now  are  gone. 
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There  is  a  consideration  of  domain,  and  we  have  a  choice  of 
procedures  in  defining  the  transformation  further.  Eqs.  (3*2)  and 
(3.9)  give 


z  =  e(A)*|  sin  k^z-Z  (3-19) 

while  also,  by  definition  of  g. 


||  =  g'  =  e(A)/|  sin  kQz.  (3-20) 

The  second  equation  shows  that  Z,  which  takes  all  values  as  8  varies 
between  any  two  successive  Integer  multiples  of  u  (cf.  eq.  (3*10)), 
rises  or  falls  monotonically  with  z,  for  fixed  A,  according  to  the 
sign  of  sin  kQa.  We  may  choose  A  >  0  for  sin  8  )>  0  and  A<  O.for 
sin  8<  0.  InUthis  case  also,  eq.  (3-19)  shows  that  z  and  Z  have  the 
same  sign  always.  There  seems  to  be  no  special  advantage  to  doing 
this,  however,  so  I  will  just  take  A  =  +  2  and  allow  Z  to  reverse 
sign  relative  to  z  each  time  8  passes  through  an  integer  multiple  of 
it.  Each  time  it  does  the  upper  (lower)  signs  in  eqs.  (3-11),  and  ^ 
(3-15)  -  ( 3 ;  1 7 )  give  way  to  the  lower  (upper)  signs.  For  finite  z 
both  Z  and  Z  are  infinite  across  these  points,  while  x  and  x  are 
continuous.  The  system  of  equations  of  the  axial  motion  along  the 
z-axis  from  -  <*>  to  +  ®,  coupled  with  the  x-motion,  has  been  replaced 
by  a  pair  of  systems  corresponding  to  motion  along  succeeding  half¬ 
cycles  of  kQz,  covering  a  combined  finite  length  27r/kQ.  And  in 
Z-space,  there  are  axially  reflecting  mirrors  at  +  ®.  The  situation 
is  illustrated  in  Figure  2. 

For  vanishing  P  ,  a  =  0  and  eqs.  (3-15)  and  (3-16)  become 
slightly  simpler,  y 

Z  -  (2C  -  x2)kQ  sinh  kQZ  =  0  (3-21) 


x  +  |  b  sech2  |  k^Z^sinh  2kgX  =  0.  (3.22) 


Note,  once  more,  the  absence  of  the  periodic  functions  of  z. 
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CHAPTER  IV 
CONCLUSION 


A  planar  undulator  wiggler  field  configuration  has  been 
constructed  and  examined  from  the  standpoint  of  a  transformation  of 
the  axial  coordinate  that  sends  the  description  of  the  single- 
electron  motion  between  z  =  +  «,  to  that  of  a  particle  in  a  box  in 
which  an  oscillation  cycle  between  Z  =  +  ®  corresponds  to  the 
traversal  of  a  single  cycle  of  wiggler  period  |Az|  =  2-n/k.Q. 

A  similar  transformation  is  possible  for  the  x-motion,  but  this 
has  not  been  discussed  as  the  resulting  system  does  not  appear 
especially  interesting. 

The  special  case  Py  *  0  should  be  an  interesting  one  physically 
(cf.  remark  in  the  paragraph  following  eq.  (2.12)),  so  eqs .  (3.21) 
and  (3.22),  in  particular,  might  repay  some  effort  of  further 
investigation. 
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